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Abstract 

Using a sample of 5.8 x 10 7 J/V> events, the Beijing Spectrometer experiment has 
searched for the decay J/ip — ► efi. Four candidates, consistent with the estimated 
background, are observed, and an upper limit on the branching fraction of J/ip — ► e/x 
of 1.1 x 10" 6 at the 90 % C.L. is obtained. 
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1 Introduction 



In the minimal standard model of electroweak theory, the three separate lepton 
numbers, electron number, muon number, and tau number, are conserved, 
but this conservation law may be broken in many extensions of the Standard 
Model, such as grand unified models [1], supersymmetric standard models [2], 
left-right symmetric models [3] and models where electroweak symmetry is 
broken dynamically [4]. Recent experimental results from Super-Kamiokande 
[5], SNO [6] and KamLAND [7] indicate strongly that neutrinos have masses 
and mix with each other; hence lepton flavor symmetries are broken. There 
have been many studies both experimentally and theoretically on searching 
for the lepton flavor violating processes [8], mainly from /x, r and Z decays [9]. 
Theoretical studies on the possibility of searching for the lepton flavor violation 
in decays of charmonium and bottomonium systems are discussed in Refs. 
[10]- [12]. In this paper we report on a search for lepton flavor violation via the 
process J/-0 -> e ± /j, T using 5.8 x 10 7 J/tp events at BEPC/BESII . 



2 BES detector 



The BEijing Spectrometer (BES) [13] is a large general purpose solenoidal 
detector at the Beijing Electron Positron Collider (BEPC). The beam energy 
of BEPC ranges from 1.0 to 2.5 GeV, and the peak luminosity at the J/ip is 
around 5 x 10 30 cm _2 s _1 . The upgrade from BESI to BESII [14] includes the 
replacement of the central drift chamber with a straw-tube vertex chamber, 
composed of 12 tracking layers arranged around a beryllium beam pipe and 
with a spatial resolution of about 90 /zm; a new barrel time-of-flight (TOF) 
counter with a time resolution of 180 ps; and a new main drift chamber (MDC), 
which has 40 tracking layers with a dE/dx resolution of &dE/dx — 8.0% and 
a momentum resolution of o v jp = 1.78% VI + p' 2 (p in GeV/c) for charged 
tracks. These upgrades augment the pre-existing calorimeter and muon track- 
ing systems. The barrel shower counter (BSC), which has a total thickness of 
12 radiation lengths and covers 80 % of Att solid angle, has an energy reso- 
lution of cte/E = 0.21/V^E (E in GeV) and a spatial resolution of 7.9 mrad 
in and 2.3 cm in z. The \i identification system consists of three double 
layers of proportional tubes interspersed in the iron flux return of the magnet. 
They provide coordinate measurements along the muon trajectories with res- 
olutions in the outermost layer of 10 cm and 12 cm in r<p and z, respectively. 
The absorber thicknesses in front of the three layers are 120, 140, and 140 
mm, and the solid angle coverage of the layers is 67 %, 67 %, and 63 % of 4n, 
respectively. 
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3 Event selection 



The initial selection of J/i/j — > e ± fi T requires that candidate events have two 
oppositely charged tracks within the polar angle region, | cos (9| < 0.8, and 
fewer than four neutral tracks. A charged track should have a good helix fit; 
a momentum, P, with 1.45 GeV/c < P < 1.65 GeV/c; and originate from the 
interaction region, defined by \x\ < 0.015 m, \y\ < 0.015 m and \z\ < 0.15 m. 
To reject cosmic rays, the TOF time difference of the two charged tracks must 
be less than 1.2 ns. The invariant mass of the e\i system, M efJi , must satisfy 
2.95 GeV/c 2 < M efI < 3.25 GeV/c 2 , and the angle between the two charged 
tracks, #i 2 , is required to be greater than 178.5°. The 9 12 distribution is shown 
in Fig. 1. 

Isolated photons are defined as those photons having an energy deposit in the 
BSC greater than 50 MeV, an angle with any charged track greater than 15°, 
and an angle between the direction defined by the first layer hit in the BSC 
and the developing direction of the cluster in the xy-plane less than 18°. There 
must be no isolated photon in the selected event. 

The above criteria select back-to-back, two prong events, such as J/t/j — > 
e + e~(7), yU + /i _ (7), K + K~ , 7r + 7r~ and e + e~ — > e + e~(7), etc. To select electrons 
and muons, the BSC and /i counter information is used. The actual selection 
criteria are based on distributions determined from data. Fig. 2 shows distri- 
butions of E/P, where E is the energy deposited in the shower counter, for 
candidate electron and muon tracks. To be an electron, a track must have no 
hits in the muon counter and satisfy E/P > 0.7. 

To select muon tracks, the differences, Si(i = x, y, z), between the closest muon 
hit and the projected MDC track in each layer are used. Fig. 3 shows these 
differences in the third layer of the \i counter. The distributions are Gaussians, 
and standard deviations, Oi{i = x,y,z), are determined for each layer of the 
fi counter. 

To reduce background, we make a tight cuts on the muon. A good hit in the 
fi counter requires \5i\ < Oi for i = x, y, and z. The total number of layers 
hit in the \x counter is denoted as /ihu and can range from zero to three. The 
number of good hits is denoted by Hf°° d - A track is considered as a muon if 
E/P < 0.3 and fx g h °° d is greater than zero if the transverse momentum of the 
track, P xy , is less than 0.75 GeV/c, pi 9 h °° d > 1 if 0.75 GeV/c < P xy < 0.95 
GeV/c, or n 9 °° d = 3 if P xy > 0.95 GeV/c. 

Using these selection criteria, four J/ip — > e// candidates are found in 5.8 x 
10 7 J/ijj events. The characteristics of these events are listed in Table 1. 
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(a) e + e"^|i + |i"(data) 
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(b) e + e"^e + e"(data) 
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(c) J/y->e|i(MC) 
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Fig. 1. The opening angle between two charged tracks for (a) e + e" — ► /U+jLi - (data), 
(b) e + e~ — > e + e _ (data) and (c) J/Y> — ► e/i (Monte Carlo data). 



4 Efficiencies and backgrounds 



The main backgrounds for J/xjj — > e/j, are nearly back-to-back J/ip —> e + e~(7), 
fi + fi~( , ~f), K + K~ , 7r + 7r _ , e + e~ — > e + e~(7), and e + e _ — > fi + fi~('~f) events where 
one or more tracks is misidentified as a /i or e. It is therefore important to 
determine the \i and e particle identification (P1D) misidentification efficien- 
cies for the background channels, as well as the PID efficiencies for the sig- 
nal channel. In this study, we determine these efficiencies using information 
from the BSC and muon counters from data. The overall efficiencies include 
the PID efficiencies and the geometric efficiencies, euc-, determined using our 
Monte Carlo (MC) program SIMBES (SIMulation at BES), which is based on 
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Fig. 2. Distribution of E/P of candidate a.) e and b.) /i tracks. E and P denote the 
energy deposited in the BSC and the momentum measured in the MDC, respectively. 



GEANT3.21. For the signal channel J/ip 
is then 



e/i, the total selection efficiency 



where e ei i-MC is the geometric efficiency, e e -). e is the electron PID efficiency, 
and e M „> M is the muon PID efficiency. Using 30000 J/ip — > e/i Monte Carlo 
events with an angular distribution of 1 + cos 2 9, the geometric efficiency is 
determined to be e eil -uc — (53.7 ± 0.3)%. 



For the background channels J/ip (e + e ) 
efficiency after J/ip — * e/i selection is 



XX, where X — e, /i, 7r, i^, the 
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Fig. 3. The differences, Sx, 5y, and Sz, between the closest hit and the extrapolated 
MDC track in the third layer of the /x system. 

where exx-MC is the Monte Carlo geometric acceptance, shown in Table 2, 
and ex->e and ex^jx are the particle misidentification efficiencies for X being 
identified as an electron or a muon, respectively. 

The e and \x PID efficiencies are determined using e-pair samples of e + e~ — > 
e + e~(7) (Bhabha, resonance, and continuum production) events and /i-pair 
samples of e + e~ — > fi + fi~('y) (resonance and continuum production) events. 
The highest momentum track is required to satisfy 1.45 < P < 1.65 GeV/c 
and to be identified as an e or /j, according to the PID selection criteria. The 
other track in the event is then assumed to be the same type of lepton, and the 
fraction of these satisfying the lepton selection determines the efficiency. The 
PID efficiencies are shown in Table 3, and the total efficiency for J/ip — > e/i 
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Table 1 



The characteristics of the J/ip — ► e/j, candidates. 9±2 is the angle between two 
charged tracks. 



RUN No. 


14676 


16419 


18940 


19149 


REC No 


16877 


26352 


27010 


9196 


/if fripVi 


S 1 1 S 

u.llO 


S 1 1 7 


S 901 


3 143 


CM 9 


179.5 


179.8 


179.6 


179.5 


track 


e fi 


e /j 


e /i 


e /i 


P (GeV/c) 


1.591 1.522 


1.570 1.543 


1.587 1.611 


1.605 1.535 


E/P 


0.7188 0.1758 


0.8577 0.1778 


0.8835 0.1415 


0.7155 0.1035 


qood 


3 


3 


3 


3 


Table 2 










Monte Carlo 


geometric efficiencies. 








channel 




MC Efficiency 






J/ip — > ee 


tee-MC 


(ct-i 47+0 09,)% 






J/V' -► MM 




(58.32±0.02)% 






J/V' — ► 7T7T 




(52.74±0.29)% 






J/V' -»• ATT 


CKK-MC 


(24.38±0.24)% 






e + e _ — > e + e~ 


(7) £ ee( 7 )-MC 


(32.51±0.03)% 






e + e~ — > /tx + /ir 


(7) e w ( 7 )-MC 


(42.96±0.29)% 





is (9.7 ±0.5)%. 



The misidentification efficiencies between electrons and muons cannot be stud- 
ied from the data because the events selected for this study are the same as 
our J/ip efi candidates. A rough estimate of this background using a Monte 
Carlo simulation predicts a total of 7 background events from e + e~ — > e + e~(7) 
and e + e~ — > n + n~{j). Since the Monte Carlo simulation is not precise enough 
to determine this background, we conservatively ignore these backgrounds in 
the estimation of the upper limit on the branching ratio for J/tp — > e/x. 

For the misidentification efficiencies of hadrons, tc and K tracks are selected 
from J/ijj — > p ± 7r T and K*^K^ decays. The track with the highest momen- 
tum is considered as the it and K in the respective samples, and a total of 
60551 7r tracks and 9275 K tracks are selected. The fraction of the tracks 
above 1.2 GeV/c passing the electron and muon selection criteria are used to 
obtain the misidentification efficiencies, that are listed in Table 3. The resul- 
tant background rates from particle misidentification of J/ip — > K + K~ , 7r + 7r~ 
events are listed in Table 4, taking into account their branching ratios from 
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PDG2002 [9]. 
Table 3 

The particle identification/misidentification efficiencies. 

regarded as e regarded as (i 

e sample 95.3%(1±0.02%) 
H sample 19.0%(1±0.6%) 
7T 3.6%(1±2.1%) 0.46%(1±5.98%) 

K 3.11%(1±5.79%) 0.38%(1±16.8%) 

Table 4 

The misidentification rates and backgrounds from hadronic channels. 

decays misidentification number of 

rate background 
J/V> -»• vttt 1.74 x 1(T 4 1.49 
J/V> 5.77 x 10~ 5 0.79 

total 2.3 

There are a number of sources of systematic error, including the error on 
the total number of J/ip (4.7%) events, the errors on the branching ratios 
from PDG2002 [9], the uncertainties of the muon and electron identification 
efficiencies, shown in Table 3, and the errors on the hadronic misidentifica- 
tion efficiencies, also shown in Table 3. Table 5 summarizes the errors in the 
background from the hadronic channels. The total background, ignoring the 
contributions from J/ip — > e + e~(7) and /i + /i~(7), is 2.3±0.3 events. Therefore 
the observed four events are consistent with background. 

Table 5 

Summary of estimated errors for hadronic background channels. 





Ntotal 


BR. 


epm 


Mpid 


total error 


number 




-> 7T7T 4.7% 


15.6% 


2.10% 


5.98% 


17.4% 


1.49±0.26 


j/y- 


+ KK 4.7% 


13.1% 


5.79% 


16.8% 


22.5% 


0.79±0.18 



5 Results and discussion 

The four observed J/ip — > e/j, candidates are consistent with the estimated 
background, hence an upper limit for J/ip — > e/x is determined. Only the back- 
grounds from the hadronic channels are taken into account in the estimation 
of the upper limit. 
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There are several ways to determine the upper limit [15] [16] [17]. The following 
formula is used to calculate the upper limit for Br(J/ip — > e/i): 

Br < \(N OB ,N BG )/[N T x e J/VMe/i ], 

where A is the upper bound of a 90% CL. for an unknown Poisson signal mean, 
for total observed events, Nob, and known mean background, Nbg- Nt is the 
total number of J ftp events, and e is the detection efficiency. \(N b, N B g) can 
be calculated using the method of Ref. [15]. 

As a conservative estimation, we take the number of background events as the 
central value of the number of background events reduced by one standard 
deviation, i.e. N B g = 2.0 and N b=^, then \(N b, N bg )=5.98. Therefore we 
obtain: 

Br(J/^j -> en) < 1.1 x KT 6 . 

In summary, we searched for J/ip — > e/i based on 5.8 x 10 7 J/ip events and 
observed four J/tp — > e/i candidates passing our selection criteria, which are 
consistent with the estimated background. The upper limit of Br(J/ip — > e/i) 
is determined to be 1.1 x 10~ 6 at the 90 % CL. 
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